Gravitational radiation, energy and reaction on quasi-spherical black holes 
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Gravitational radiation is locally defined where the wave- 
fronts are roughly spherical. A local energy tensor is defined 
for the gravitational radiation. Including this energy tensor as 
a source in the truncated Einstein equations describes gravi- 
tational radiation reaction, such as back-reaction on a roughly 
spherical black hole. The energy-momentum in a canonical 
frame is covariantly conserved. The strain to be measured by 
a distant detector is simply defined. 
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Gravitational-wave astronomy is expected to become 
a major observational science in the coming millcmum. 
Gravitational-wave theory, however, is based on various 
approximations, as there is no general local definition of 
a gravitational wave in General Relativity; see e.g the 
review of Thorne [jy. Einstein Q initiated the study of 
gravitational waves, in the linearized approximation on 
flat space-time, but there were objections to their reality 
and ability to carry energy. General belief came only in 
the 1960s with two further approximations. Firstly, as- 
suming asymptotic flatness allowed a definition of gravi- 
tational wave infinitely far from all sources, as initiated 
by Bondi and formalized in the conformal framework 
of Penrose J3|. Secondly, the high-frequency or short- 
wavelength approximation allowed a locally averaged def- 
inition of gravitational wave, as formulated by Isaacson 
after initial ideas of Wheeler In both cases, the 

physical reality of the waves is indicated by definitions of 
their energy which, added to that of the matter, yields 
energy-balance equations: the Bondi energy- loss equa- 
tion and the averaged Einstein equation, respectively. 
This letter describes how this may also be achieved in 
a quasi-spherical approximation scheme. Moreover, the 
resulting energy tensor of the radiation is local. Thus the 
approximation allows a definition of gravitational radia- 
tion at any point. 

The recent quasi-spherical approximation || was in- 
tended to be applicable to coalescing black-hole bina- 
ries, one of the main expected sources for gravitational- 
wave detectors. It may also prove applicable to neutron 
stars or supernovas, or indeed any astrophysical situation 
which has rough spherical symmetry. Further, for any 
other process which can be enclosed by roughly spherical 
surfaces, it provides a mid-zone and, assuming isolation, 
far-zone approximation. Mathematically this is achieved 



by linearizing certain fields which would vanish in ex- 
act spherical symmetry, having made a decomposition 
of space-time adapted to the roughly spherical surfaces, 
henceforth called transverse surfaces. Since this makes 
no assumption of closeness to stationarity, the approx- 
imation holds for arbitrarily fast dynamical processes. 
The approximation has also been tested against angular 
momentum by applying it to Kerr black holes : the er- 
ror in the strain waveform is much lower than expected 
signals from binary black-hole coalescence. 

The wavefronts of outgoing and ingoing gravitational 
radiation form two families of null hypersurfaces, inter- 
secting in the two-parameter family of transverse spatial 
surfaces. This geometry is described by the formalism of 
dual-null dynamics |^|,^), briefly summarized here. La- 
belling the hypersurfaces by null coordinates x + and x~ , 
and taking coordinates x a for the transverse surfaces, the 
space-time metric takes the form 



g = h ab (dx a + s a + dx + + s a _dx ) <g) 

(dx b + s\dx + + s h _dx~) - 2e~ f dx + <g> dx~ 



(1) 



where h is the metric of the transverse surfaces, s± are 
two shift vectors and / is a normalization function; where 
appropriate, Latin letters are used for transverse indices 
and Greek letters for space-time indices. The covariant 
derivative of h is denoted by D, and the evolution deriva- 
tives are 



A+ = ±L 



d/dxl 



(2) 



where L denotes the Lie derivative and _L denotes pro- 
jection by h. It is also useful to decompose h into a 
conformal factor Q and a conformal metric k by 



(3) 



such that the Hodge operator * of k is given by the stan- 
dard spherical-polar area form 



il = sin 6 d9 A d<p 



(4) 



where x a = (0, 4>) are quasi-spherical polar coordinates. 
Then f2 _1 is the quasi-spherical radius. The extrinsic 
curvature of the dual-null foliation may be encoded in 
conformally rescaled expansions §± and traceless shears 
<r± a b, inaffinitics v± and twist uj a , defined previously 0] 
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and below as necessary. The Einstein equation and con- 
tracted Bianchi identity may then be written in a first- 
order form expressing A± derivatives of the dynamical 
fields v±, tt, lj, Cl, f, k, s±) in terms of the dynami- 
cal fields and their first and second transverse derivatives. 

The first quasi- spherical approximation, introduced 
previously |J, consists of linearizing in 0;±,u),8±,D), 
which vanish in spherical symmetry [pL0|Jl if] . This yields 
a greatly simplified set of truncated equations, decou- 
pling into a three-level hierarchy. Moreover, the last 
level, the equations for (ui,s±), need not be solved for 
the radiation problem. This is because, fixing A + to be 
the outgoing derivative, the Bondi news at null infinity 
3 ± is essentially ? T = Q,~ 1 A T k [|6|j7|Jl^|. (Generally, one 
should use the null derivatives A± — -LL S± ; however, ap- 
plied to transverse tensors, they coincide with A± in the 
quasi-spherical truncation, so will not be distinguished in 
this article). Specifically, the conformal strain tensor at 
future null infinity S + is 



T7 / <!■'■ 



This means that 



£ 

e = - 

r 



(5) 



(0) 



is the strain tensor at a large distance r from the 
source, so that the displacements to be measured by a 
gravitational-wave detector are 



S(: 



e ab e a e» 



(J) 



where the Cartesian basis vector e is the direction of dis- 
placement. Thus the variables of the quasi-spherical ap- 
proximation are directly related to the observable strain; 
no further far-zone approximation is required. 

Moreover, the Bondi flux may be localized in terms 
of a gravitational-radiation energy flux, derived below as 
the appropriate contraction of an effective energy tensor 
of the gravitational radiation: 







a/3 



(A a k,A f3 k) - yi s {A 7 k,A 5 k)g af , 



32tt 



(8) 



where (a, (3) = k ab k cd a ac pbd, A a = and units are 
such that Newton's gravitational constant is unity. Ex- 
plaining this is the main purpose of this letter. Princi- 
pally, acts like a matter energy tensor in the truncated 
Einstein equations and satisfies a covariant conservation 
law, conservation of energy. Also, O satisfies the strong, 
dominant, weak and null energy conditions fl3|| . Thus 
gravitational radiation carries positive energy. It should 
be stressed that (i) mathematically, is a genuine ten- 
sor, but depends on the dual-null foliation, not just on 
the space-time; (ii) the physical interpretation of as 
energy requires the quasi-spherical approximation to be 
valid, meaning that the transverse surfaces must indeed 



be roughly spherical. This will not be made precise here, 
as the range of validity of the approximation is not clear 
in advance and best explored in applications. The intu- 
itive meaning of roughly spherical should be clear by any 
standards. 

In short, the quasi-spherical approximation allows a 
local definition of the energy-momentum-stress of grav- 
itational radiation, and therefore of the radiation itself: 
gravitational radiation is present at a given point if and 
only if is non-zero there. With the above orientation, 
there is outgoing radiation if and only if A_fc (equiva- 
lently is non-zero, and ingoing radiation if and only 
if A + fc (equivalcntly is non-zero. The terminology 
radiation instead of wave is generally preferable, since 
A±k need not be oscillatory. Instead, frequency spectra 
may be defined by Fourier transformations of k from x 
(at constant x^) to frequency f± [fL3f . Then a gravita- 
tional wave is gravitational radiation which is peaked in 
frequency space. This reflects the different physical ba- 
sis as compared to that of the Isaacson high-frequency 
approximation |,§. In the quasi-spherical approxima- 
tion, gravitational radiation is defined even if there is no 
typical frequency. 

Writing the non-zero components 



6±± = ||A ± fe||7327r 

e ab = e^(A + k 1 A_k)h ab /32TT 



(9) 
(10) 



where ||a|| 2 = {a, a), the 0±± components are remi- 
niscent of the high-frequency approximation. The Q ab 
components show that there is generally a transverse ra- 
diation pressure produced by a combination of ingoing 
and outgoing gravitational radiation, as well as the ex- 
pected radial radiation pressure. The fact that 0-i = 

may be interpreted as meaning that gravitational radia- 
tion is purely radiative and workless, as for the massless 
Klein-Gordon field. 

Apart from the fact that a gravitational wave generally 
has two polarizations, as encoded in the two independent 
components of k, is analogous to the energy tensor of 
Einstein- Rosen gravitational waves (lj] , which takes the 
massless Klein-Gordon form in terms of a gravitational 
potential generalizing the Newtonian potential. A com- 
plex gravitational-radiation potential may similarly be 
defined for 0, which then takes a scalar-field form p3| . 

If is to measure the energy of the gravitational ra- 
diation, then gravitational radiation reaction, the back- 
reaction on the space-time, should be described by in- 
cluding as a matter energy tensor in the truncated 
Einstein equations. In fact, there is a more logical way to 
formulate this, suggested by noting that the gravitational 
radiation is encoded in the conformal shears <r±, which 
are linearized in the first approximation. Then retaining 
non-linear terms in c± should give a more accurate ap- 
proximation for the gravitational-radiation sector of the 
theory. Thus the second quasi- spherical approximation, 
described in detail in a longer article [jl3| , consists of lin- 
earizing in (u>, s±,D) only. The resulting equations also 
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decouple, this time into only two levels, with the last level 
for {lo,s±) again being irrelevant to the radiation prob- 
lem. The remaining equations are, taking the vacuum 
case, 



ft 2 i? 



± 



A ± n= .,. 
A±k = 

A±0± = -i/±tf±-±fi|k±|| 2 

A±tf T = + e _/ ) 

A±u T = -n 2 (^+tf- + e-f - *(?+,?_>) 



where (a o /3) a f, = fc cd ' 



(11) 
(12) 

(13) 
(14) 
(15) 
(16) 
(17) 



"Oiacflbd- The first three equations 
effectively define ($±, v±, <;±). The first approximation 
may, of course, be recovered by linearizing in <;±. The 
additional quadratic terms in q± appear in the same way 
that a matter energy tensor with the form of O would, 
as can be seen by comparing with the spherically sym- 
metric equations @|J. More formally, one may intro- 
duce a truncated Einstein tensor C |13| in the first ap- 
proximation. Then the truncated Einstein equations are 
C = 8nT in the first approximation, where T is the en- 
ergy tensor of the matter, and C = 87r(T + 9) in the 
second approximation. This is the first reason for iden- 
tifying as an energy tensor for the gravitational radia- 
tion: it plays the role of an effective matter energy tensor 
in the second approximation. 

Mathematically, the difference between first and sec- 
ond approximations is that in the first approximation, 
the equations for (i?±, v±, f2, /), the variables which sur- 
vive in spherical symmetry, decouple from the equations 
for (?-t, k), which constitute a wave equation for k. Phys- 
ically this describes gravitational-wave propagation on a 
quasi-spherical background. The background is not fixed 
in advance and need not be spherically symmetric, so 
even the first approximation is widely applicable. There 
is no such decoupling in the second approximation: the 
gravitational-radiation terms (<^±, k) now enter the equa- 
tions for the quasi-spherical part of the geometry, which 
thus reacts to the passage of the waves. Gravitational 
radiation reaction has thereby been included; there is no 
longer a background which is independent of the waves. 

Nevertheless, both first and second approximations 
share the remarkable feature that, to compute the observ- 
able waveforms, no transverse derivatives need be consid- 
ered. The truncated equations form an effectively two- 
dimensional system, to be integrated independently at 
each angle of the sphere. Physically this means that the 
observed gravitational-wave signal depends only on the 
line of sight to the source, surely a plausible result. More- 
over, by virtue of the dual-null formulation, the equations 
are already written in characteristic form, the mathemat- 
ically standard form for analysis of hyperbolic equations. 
Numerical implementation is consequently straightfor- 
ward and computationally inexpensive. Numerical codes 
exist for both first and second approximations 0. 



The second reason for identifying 9 as an energy tensor 
is that, added to the energy tensor of the matter, it yields 
a covariant energy conservation law. This is derived in 
general in the longer article fis}] and in the vacuum case 
as follows, by a more direct method. As in spherical 
symmetry ]TqJii| , there is a canonical flow of time defined 
by the vector or 1-form 

£ = *dfT 1 



(18) 



where * is the Hodge operator of the evolution space, 
*1 = e~?dx + A dx~ . Its non-zero components are 



(19) 



Then £ is analogous to the Killing vector of a stationary 
space-time. The energy-momentum density of the grav- 
itational radiation, referred to the canonical flow, is the 
vector or 1-form 



with non-zero components 



J± 



= ±e f e±±A T n 



(20) 



(21) 



This has the same form as the energy-momentum density 
of the matter in spherical symmetry pO 11 1 . The corre- 
sponding energy flux of the gravitational radiation is the 
dual 1-form 



t/> = *j. 

The conformally rescaled flux 

ip = fr 2 ?A 

then has non-zero components 



<p± 



=fIIg=I 



64tt 



(22) 



(23) 



(24) 



These expressions have the same form as those for the 
Bondi flux at 3 T pj7|,p^|. Thus the conformal flux <p is 
a local generalization of the Bondi flux. Denoting the 
Hodge operator of g by *, the quasi-spherical truncation 
identity 



*i = *i a *rr 2 

allows the divergence of j to be written as 

V Q j" = -kdkj = f7 2 *d(fT 2 *j) = tt 2 *dtp. 
The truncated Einstein equations dO-llJ) yield 



ilk+IHMI 2 ) 



and therefore 

*dtp = e^(A_tp + — A + <yS_) = 0. 



(25) 



(26) 



(27) 



(28) 
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Thus j is covariantly conserved: 

v a j a = o. 



(29) 



Physically this represents conservation of energy, as in 
spherical symmetry jl0],[ll]]. In the presence of matter, 
it holds for the combined energy of the gravitational ra- 
diation and matter |^3[ , as in cylindrical symmetry [p"4[ . 
The Noether charge associated with the Noether current 
j provides a definition of active gravitational mass-energy 
generalizing that of spherical symmetry, including the en- 
ergy of the gravitational radiation. Then energy conser- 
vation can be written in the form of a first law |l3| ] . 

In summary, the quasi-spherical approximation scheme 
allows local definitions of gravitational radiation and its 
energy. The radiation is superimposed on a background 
in the first approximation, but reacts back on the space- 
time in the second approximation. Comparing results of 
the first and second approximations, in a given situation, 
provides a valuable internal guide to their accuracy, inde- 
pendent of any other estimates. In particular, the results 
obtained for Kerr black holes are numerically indistin- 
guishable in the first and second approximations . 

Gravitational-wave theory has progressed since the re- 
view of Thome M by higher-order post-Newtonian ap- 
proximations |p.5| and close- limit approximations [pL6|| , 
which can be used to describe, respectively, the pre- 
coalescence and post-coalescence phases of a binary 
black-hole inspiral. The coalescence phase has long been 
thought to be tractable only by numerical methods, but 
despite great efforts [|L7J, the desired waveforms are not 
yet known. The quasi-spherical approximation is in- 
tended to apply to the post-coalescence phase, thereby 
reducing the period for which full numerical simulations 
are required. Similarly, a quasi-equilibrium approxima- 
tion has also recently been suggested for the later pre- 
coalescence phase Q . 

Gravitational radiation reaction is an essential physical 
ingredient in the post-Newtonian and quasi-equilibrium 
approximations, but is not included in the close-limit ap- 
proximation, or indeed any other perturbative approach 
with a fixed background. In comparison, the second 
quasi-spherical approximation allows a fully relativistic 
inclusion of radiation reaction for rapidly evolving black 
holes. It is applicable from the very moment of coales- 
cence, meaning the appearance of trapped surfaces en- 
closing both original trapped regions. 

As an illustration, the following scenario is plausible 
on physical grounds: a distorted black hole generally 
emits gravitational radiation, which generally backscat- 
ters to produce ingoing radiation, which is absorbed by 
the black hole, thereby increasing its mass and area and 
generally changing its shape, thereby also absorbing some 
outgoing radiation, thereby changing the profile of subse- 
quent gravitational-radiation emission, and so on by feed- 
back. This entire process is actually described, step by 
step, by the truncated Einstein equations (p~Tj-|l7|) . The 
area-increase property was derived as the second law of 



black-hole dynamics JlSfl , where a black hole is locally 
defined by a type of trapping horizon, which was also 
shown to be achronal and have topologically spherical 
sections. Corresponding spherically symmetric first |Tl| l 
and zeroth p(| laws of black-hole dynamics admit quasi- 
spherical generalizations (l3| involving a local definition 
of surface gravity and the above mass. 

Thus there is an astrophysically realistic context in 
which both black holes and gravitational radiation are 
locally defined, along with their physical attributes, with 
each influencing the other. Apart from practical appli- 
cations to gravitational-wave astronomy, this provides a 
rich arena in which to advance theoretical understand- 
ing of the dynamical interaction between gravitational 
radiation and black holes. 
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